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We consider soliton creation when there are "twist" degrees of freedom present in the model in 
addition to those that make up the soliton. Specifically we consider a deformed 0(3) sigma model 
in 1+1 dimensions, which reduces to the sine-Gordon model in the zero twist sector. We study the 
scattering of two or more breather solutions as a function of twist, and find soliton creation for a 
range of parameters. We speculate on the application of these ideas, in particular on the possible role 
of magnetic helicity, to the production of magnetic monopoles, and the violation of baryon number 
in nuclear scattering experiments. 

Field theories are known to have two sectors - small excitations called "particles", and large excitations called 
"solitons" . Perturbative quantum field theory deals with particles and their interactions. Solitons are generally studied 
as classical solutions, unrelated to the particle description. An outstanding open problem in quantum field theory is 
to describe solitons as an assembly of particles. A closely related problem is the central focus of this paper - how can 
we build solitons by scattering particles? 

The relation between solitons and particles is explicitly known in the sine-Gordon model [1] where the soliton 
operator is given in terms of an infinite number of particle operators. Thus to build solitons from particles in the sine- 
Gordon model requires assembly of an infinite number of particles and is impossible. This ties in with the complete 
integrability of the sine-Gordon model wherein soliton number is preserved. Therefore it is necessary to consider 
models that are not completely integrable. However integrable models are amenable to a lot of analysis, and it seems 
unwise to abandon them altogether. We would like to take advantage of the insights that are provided by integrable 
models and to use them in a non-integrable setting. Hence we aim to start with a completely integrable model, modify 
it suitably to eliminate complete integrability while still preserving as many features as possible, and then study 
soliton creation in the non-integrable model. 

There can of course be many modifications of a completely integrable model. Indeed, the A^ 4 model in 1+1 
dimensions studied in [2-4] can also be viewed as a deformation of the sine-Gordon model. Here we would like to 
study a different kind of modification; one that may be useful in generalization to higher spatial dimensions where 
we would like to study the creation of vortices and magnetic monopoles. The modification is inspired by recent work 
within the electroweak model [5, 6]: the annihilation of "twisted" magnetic monopoles leads to the production of 
magnetic fields that inherit the twist in the form of magnetic helicity. We will describe this observation in greater 
detail in Sec. IV. For the time being we note that magnetic monopoles can carry a relative twist, and so we would like 
to introduce a similar twist degree of freedom in the 1+1 dimensional sine-Gordon model, thus also breaking complete 
integrability. Then the kinks of the sine-Gordon model will be able to carry relative twist, and solitons may be created 
in scattering that involves the twist degree of freedom. The hope is that twisted initial conditions are better suited 
to the creation of solitons. In our toy model in Sec. I we will see that twist is essential for the creation of solitons 
because untwisted initial conditions lie in a completely integrable sub-space of the model. 

A different view of the soliton creation problem also suggests that twist can play a useful role. When two particles 
scatter, it is advantageous if they spend a long time together, so that other particles also have time to interact and 
create a multi-particle state such as the soliton-antisoliton pair. In our case, we will take the initial state to be a 
sequence of sine-Gordon breathers - the closest we can come to a sequence of particles in classical theory. When two 
in-phase breathers scatter in the sine-Gordon model, they pass through each other but with a negative time delay i.e. 
a time advance. This is because of an attractive force between the breathers. If the breathers are twisted, however, 
the force can be repulsive, and the scattering can lead to a time delay. In Sec. II we will study the scattering of two 
breathers as a function of twist. The numerical results confirm that twist can lead to a positive time delay in the 
scattering 

The paper is organized as follows. In Sec. I we describe our toy model, and various solutions in it. In Sec. II we 
study the scattering of two breather solutions of the model and show that twisted breathers lead to a time delay in 
the scattering. As a special case we find the time delay in the scattering of two sine-Gordon breathers. In Sec. Ill we 
study the scattering of two sequences ( "trains" ) of twisted breathers by numerically evolving the equations of motion. 
Here we find examples of soliton creation. We discuss generalizations to higher spacetime dimensions, and to vortices 
and monopoles in Sec. IV. The production of magnetic monopoles is intimately related to baryon number violation 
via electroweak sphaleron production. We discuss this connection in Sec. IV, where we also conclude. In Appendix A 
we provide some notes related to the numerical evolution of the constrained 0(3) system. 
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I. DEFORMED 0(3) MODEL 

With the motivation described above, we now move on to describe a model in 1+1 dimensions that contains kinks 
and also allows for twist. We start with the 0(3) sigma model given by the action 

S, = I #x\{d»hf (1) 



2 

where n = h(t, x) satisfies the constraint 

n 2 = 1 (2) 
This model is known to be completely integrable [7] and we will need to spoil this feature. So we consider 



S = / dfx 



(3) 



where n 3 is the third (or z— ) component of h. The additional potential term spoils the 0(3) invariance and we now 
only have Zi symmetry under 713 — > — 713 and 0(2) symmetry under rotations about the z-axis. The true vacua are 

h = ±(0,0,1) (4) 

The model has kink solutions that interpolate between the two vacua and there is a one parameter family of such 
kinks. 

One way to proceed is by introducing a Lagrange multiplier, £, to enforce the constraint n 2 = 1 



L=\{d,nf~^-{l-n 2 )- l -[l „i C, 

This leads to the equation of motion 



□n + {d^n) 2 fi — n 3 (e3 — n 3 n) = 0, n 2 = 1 (6) 
To gain some intuition, we consider an explicit solution to the constraint equation 

n = (sin 6 cos (j>, sin 9 sin cj), cos 9) (7) 



Then 



and the equations of motion are 



where a = 19. 

If <fi is constant, then 



L = l(d^9) 2 + i sin 2 9{8^) 2 - t sin 2 9 (S) 



□a + sina(l- (d^) 2 ) = (9) 
M (ffln 2 0fl"0) = (10) 



□a = — sin a (11) 



which is the equation of motion for the sine-Gordon model, and all solutions of the sine-Gordon model are also 
solutions of our deformed 0(3) model. Hence our model contains embedded spaces given by <fi = constant where the 
dynamics is that of the sine-Gordon model 1 . Thus our model contains a one parameter family of sine-Gordon models, 
one per circle of longitude; the other degree of freedom, </>, will be referred to as the "twist" . 
The kink solution in our model is exactly the same as in the sine-Gordon model 

9 k = 2tan~ 1 [e a; ], = constant (12) 



1 In quantum theory fluctuations in the <f> direction will contribute to the dynamics of 6. 
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The energy in the fields is given by 

dxH(t,x) = dx 1 h sin 2 6— + sin 2 9- h - sin 2 I 

J — 2. 2 2 2 

and for the kink evaluates to 



(13) 



E k = 2 (14) 

Boosted breather solutions can also be borrowed from the sine-Gordon model 

•q sin(wT) 



9b(t, x; xq, v) = 2 tan 1 



= constant (15) 



cosh^wX) 

where w G [0, 1] is a parameter - the oscillation frequency of the breather - and 

ij = w-Vl-w 2 • (16) 

The boosted coordinates (T, X) are given by 

T = 7 [t - v(x - x )], X = j[(x - x ) - vt] (17) 



where xq is the location of the breather at t = 0, v is the speed of the breather, and 7 = — v 2 is the Lorentz 

boost factor. 

The energy of a boosted breather is 

E b = 47^1 -uj 2 (18) 

Additional solutions can be obtained if 9 M = fc M , a constant two vector, then again 

\3 y a = — sin a (19) 

where the D'Alembertian operator now contains derivatives with respect to ?/' = \/l — k 2 x^. The <j) equation is 
satisfied provided k^d^O — 0. This can be arranged by choosing 9 = 9(p fl x fl ) with k ■ p = 0. For example, if we take 
k^ = (k,0), the solution is "spinning" in field space. In this case, 9 can be any static solution of the sine-Gordon 
model, in particular, it can be the sine-Gordon kink. This solution can be pictured on the two-sphere as lying on a 
line of longitude that spins with angular frequency k. 

The feature that our model has sine-Gordon subspaces is advantageous becase breathers are non-dissipative solutions 
in the model. Thus we can formulate non-dissipative initial conditions by writing them in terms of a sequence of 
breathers. In the quantized sine-Gordon model, the smallest breather corresponds to a particle, and hence, a sequence 
of classical breathers is the closest we can come classically to a sequence of quantum particles [8] . 

II. SCATTERING OF TWO BREATHERS 

We first study the scattering of two breathers and find the time delay in the scattering. In the untwisted case, 
this is well-defined because the breathers pass through each other. In the twisted case, the breathers get deformed on 
interaction. However, we can still calculate the time delay in the propagation of the center of energy. We will elaborate 
on this below. 

The initial conditions corresponding to two incoming breathers with twist £ are 

9(t = 0,x) = 9 b (t = 0,x;-x ,+v) + 6 b (t = 0,x;+x ,-v) (20) 

(/>(t = 0,x) = tanh(x/w) (21) 

9(t = 0,x) = 9 b (t = 0,x;-x ,v)+9 b (t = 0,x;+x ,-v) (22) 

0(t = O,x) = (23) 

where 9 b (t — 0,x;xq,v) is Eq. (15) evaluated at t = and 9 b (t = 0, x;x ,v) is obtained by differentiating Eq. (15) 
with respect to time and then evaluating at t = 0. The expressions are sufficiently messy that we don't display them. 
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FIG. 1: Time delay vs. twist for uj — 0.86 and v ranging from 0.1 (dashed curve with larger amplitude) to 0.9 (dashed curve 
with smaller amplitude). The time delay for v = 0.9 is at the 10~ 2 level though the dependence on twist is similar to that of 
the other curves. 



The <j) field is taken to have a tanh profile with a width w that is much smaller than xq. Since 9(0, x) ~ near x = 0, 
the precise choice of w is not significant. 

One subtlety arises since, in spherical coordinates, 9 £ [0, 7r] and 4> £ [0, 2ir], whereas in Eq. (20) 9 can become 
negative. This will not be of consequence, however, since we will always work with the vector n using Eq. (7), and 
dffi. In the construction of h and dth, negative 9 yields the same vector as positive 9 but with <p —></> + ir. 

The domain of the twist parameter £ is determined by noting that the total change in <f> is 27r£. Hence it is sufficient 
to take £ € [0, 0.5], since £ = 0.5 corresponds to 4>{x = — oo) = — 7r/2 and <j>(x = +oo) = +ir/2 i.e. the dynamics occurs 
on the great circle in the yz— plane, where it is given by the sine-Gordon model. For £ = 0, the incoming breathers 
are in phase; for £ = 0.5, the incoming breathers are out of phase. 

The numerical algorithm is described in Appendix A. To find the time delay, we evaluate the position of the center 
of energy on half the space 

*CE(t) = -i / dxxH{t.x) (24) 
h / 2 Jo 

where H is the Hamiltonian density (Eq. (13)) and E is the total energy. The time delay after some long time T is 
given by 

5t = T _XCE(T) (25) 
V 

This formula assumes that all the energy in the final state is propagating with speed v and so the time delay makes 
strict sense only if the incoming breathers pass through each other after interacting and retain their identity. This 
will only happen for £ = 0, 0.5 when the dynamics is sine-Gordon. For other values of £, the scattering causes some 
radiation, though most of the energy still resides in a breather-like lump and the definition in Eq. (25) can be used 
as a measure of the time delay. 

In Fig. 1 we plot the time delay versus twist for fixed uj = 0.86 and for incoming velocity ranging from 0.1 to 0.9. 
The plot shows that there is an equal time advance for the in-phase (£ = 0) and out-of-phase (£ = 0.5) sine-Gordon 
breathers 2 . The plot also shows that twist can cause a time delay and the effect is largest at low velocities. Similarly 



2 Here we disagree with Ref. [9] where a time delay is found for out-of-phase scattering on the basis of a proposed collective coordinate 
method. 
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FIG. 2: Time delay vs. twist for v — 0.1 and ui ranging from 0.86 (dashed curve with larger amplitude) to 0.95 (dashed curve 
with smaller amplitude). 



in Fig. 2 we see the role of twist on the time delay for a variety of breathers of different frequency, w, but with fixed 
v = 0.1. The time delay is largest for the low frequency (i.e. large amplitude) breathers. 



III. SCATTERING OF MANY BREATHERS 



Now we consider the scattering of two trains of N identical breathers. The total incoming energy has to be larger 
than the kink-antikink energy. With the kink and breather energies in Eqs. (14) and (18), we get the condition 

8iV 7 \/l - w 2 > 4 (26) 

which translates into 



1 l-v 2 

This condition is simply an energy requirement. 

The initial condition is a generalization of Eqs. (20)-(23) 

(28) 
(29) 
(30) 

(31) 
with 

Xj = x + (j - I) a (32) 

where Xq is the location of the innermost breather and a is the spacing between the breathers in a train. An example 
of the initial condition is shown in Fig. 3. 







N 




6(t = 


Q,x) 


= ^ [6 b (t — 0, x; —xj, +v) - 


- 6f,(t = 0, x; +Xj, —v 


<Kt = 


0,x) 


— 7r^ tanh(x/w) 

N 




6{t = 


0,x) 


= [^C* = 0,x;-Xj,v) + 

3 = 1 


9 b (t = 0, x; +Xj,—v) 


4(t = 


0,x) 


= 





6 





v 1 y "7" 


IWWW 




wwwi 




0.5 










- 


1 

e 










_ 


-0.5 












-1 


-200 





200 



X 



FIG. 3: Three snapshots of the z— component of h for the run with v — 0.5, f = 0.1 and the remaining parameters as in 
Eq. (33). The initial breather state is the black solid curve. The final state is the solid (red) curve with the kink and the 
antikink interpolating between 713 = +1 and 713 = — 1. An intermediate state is also shown (dashed blue curve). 

With the initial conditions in terms of 9 and we construct the vector field n(0, x) and its time derivative, dth(0, x). 
We then evolve h using Eq. (6), as discussed in greater detail in Appendix A. We hold the following parameters fixed 
in our numerical runs 

4 2 

AT = 4, w = 0.95, x = 32, a= = 11.1, w= = 5.55, (33) 

and v € [0, 1] and £ € [0, 0.5] are scanned over. Most runs do not yield soliton creation but some runs do. An example 
of soliton creation occurs for v = 0.5 and £ = 0.1, and is shown in Fig. 3. In Fig. 4 we mark the points in parameter 
space where we have found soliton creation. Note that non-zero twist is essential for soliton creation in this model. 

It is worth pointing out that the outcome is sensitive to the choice of our fixed parameters in Eq. (33). Dependence 
on N and the spacing a is expected; but we also find dependence on xq , suggesting that the phase of the breather trains 
on collision may be important in the outcome. This would be consistent with the chaotic nature of kink production 
in X<f> 4 observed in [2, 3]. In future work we are planning a more extensive scan of parameter space and to examine 
the chaotic nature of the scattering more closely [10]. 

IV. GENERALIZATIONS AND CONCLUSIONS 

The deformed 0(3) sigma model suggests many generalizations that all seem worthwhile to explore. 

The first generalization is to stay in 1+1 dimensions and to extend the model to an N— component vector. (This 
N is not to be confused with the number of breathers in the previous section.) It may be possible to make some 
analytic advances by considering the kink creation problem in the large N approximation. Other generalizations are 
to complexify the vector n and consider the CP(N) model, and to also include gauge fields. 

The sigma model framework can be used to generalize to vortex and monopolc formation, by simply choosing 
different potentials. We have already studied kink creation in the 0(3) model with the potential 

V{ti) = \[l-nl] (34) 
To study vortex formation in 2+1 dimensions, it would be natural to choose 



V(n) = -[l-(n 2 1+ nl)} 



(35) 
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FIG. 4: Diamonds indicate points in the velocity-twist plane where we have found soliton creation to occur with parameters 
as in Eq. (33). The point denoted by the filled circle corresponds to the parameters used in the successful kink creation shown 
in Fig. 3. Our numerical strategy scans in parameter space in steps of 0.05 in twist in the interval [0,0.5], and steps of 0.05 in 
velocity in the interval [0.4,0.95]. A more refined scan would fill in some of the parameter space though we expect a fractal 
distribution of parameters which successfully produce kinks [2, 3] 



Now the true vacua lie around the equator of an S" 2 , and the vortex solutions correspond to either the northern or 
the southern hemisphere, and the center of the vortex will be either at the north or south poles. These two solutions 
are related by the Zi symmetry under n% — > —n^. In this model, we do not have a "twist" degree of freedom as in the 
kink case. Recall that in the kink case, we had a 1-parameter family of kink solutions labeled by the constant value 
of the azimuthal field cj> as in Eq. (12). Similarly we need a one parameter family of vortex solutions and this can be 
obtained by extending the model to 0(4). Then the vacua are still given by circles n\ + n?, = 1 but there is a one 
parameter family of possible vortex solutions labeled by the angle in the (713,^4) plane. 

A point to remember about vortices is that the energy of a global vortex-antivortex pair grows logarithmically with 
separation. So their creation is necessarily transient, unless the model is gauged. A similar comment applies to the 
creation of magnetic monopoles. 

To study monopole formation in 3+1 dimensions in our 0(3) model we would choose the potential, 

V{h) = \[l-{nl + ni+nl)\ (36) 

except that this vanishes due to the constraint. The 0(3) model does have the correct topology (112 (S 2 ) = Z) for 
magnetic monopoles, but the monopoles are singular because there is no escape from the vacuum manifold. If we 
choose the potential in Eq. (36) but in an 0(4) model, we get a situation similar to the vortex discussed above in 
the 0(3) model: there are two non-singular monopole solutions corresponding to the two "hemispheres" in S 3 . To 
include a twist similar to the one for kinks, we should go to at least an 0(5) model, where the one parameter family 
of monopole solutions is labeled by the angle in the (71.4, 71-5) plane. However, monopoles can be relatively twisted even 
in the 0(3) model and we now discuss this feature in greater detail. 

Let us assume that a monopole is located at z = +a and its asymptotic vector field n is given by Eq. (7) with 
8 = 6+, where 8+ is the spherical angle measured from the +z— axis but with the origin of the coordinate system at 
z = +a. An antimonopole is now placed at z — —a and its asymptotic field is given by 

n = (sin 8 cos <f>, sin 8 sin (f>, cos 8) (37) 

with 8 = 7T — where is the spherical angle measured from the +z— axis but with origin at z = —a. The azimuthal 
angle (f> = <p + a and a plays the role of the twist in the 0(3) model. 

In Rcf. [11] Taubes used a Morse theory construction in the gauged 0(3) model and proved the existence of a static 
solution consisting of a monopole and an antimonopole (mm) separated by some distance. This is surprising because 



one expects the Coulombic — 1/r 2 attractive force between the mm) to bring them together, after which the two 
can annihilate. In Taubes' solution, the mm are kept apart by a repulsive force provided by the "twist" . The result 
is an unstable solution since the monopoles can untwist, come together, and annihilate. The solution itself is quite 
complicated but has been found numerically in Refs. [12, 13]. 

A construction similar to the one used by Taubes exists in the electroweak model and the resulting mm solution is 
called the "sphaleron" [14, 15] (for a review see [16]). The form of the sphaleron solution is spherically symmetric for 
a particular value of coupling constants and is known explicitly up to two simple profile functions given by ordinary 
second order differential equations. The solution has been studied numerically for a large range of coupling constants 
[17-19]. 

A particularly convenient feature of the Morse theory construction, is that it provides a path from the mm solution 
to the vacuum. Thus the sphaleron decay path describes the annihilation of a twisted mm pair. As shown in Refs. [5, 6] 
the decay products of the annihilation include electromagnetic fields with non-zero magnetic helicity 3 , which is defined 
as 



H = 



d 3 xA ■ B 



(38) 



The magnetic helicity, which measures the twisting and linking of magnetic field lines, is inherited from the initial 
twist in the mm pair that make up the sphaleron. The sphaleron decay analysis also shows that, after a short period 
of time, the evolution of the electroweak fields is such that magnetic helicity is conserved. 

By a time reversal of the sphaleron decay process, we conclude that the convergence of electromagnetic helicity 
can possibly lead to the creation of sphalerons, and hence also mm with a relative twist. The twist may also provide 
the necessary repulsive force between monopoles and antimonoples so as to separate them and to prevent them from 
re- annihilating. 

The usefulness of magnetic helicity may be conjectured based on another line of reasoning. In highly conducting 
plasmas it is known that magnetic helicity is a conserved quantity. The derivation goes as follows: 



<m 

dt 



J d 3 x [ d t A • B + A ■ d t B ] 
J d 3 x [(E + VA°) • B - A ■ V x E] 
f d 3 x [2E • B + (VA°) • B] 
J d 3 x [2E • B - ,4°V ■ B] 



2 / d 3 a 



v x B • B] 



-> 



(39) 



In the last two steps we have used Ohm's law as applied to a plasma element moving with velocity v. The electrical 
conductivity is assumed to be infinite, a — » oo, though there are arguments in support of helicity conservation even 
if a is finite but large. Crucially for us, in the second and fifth steps, we have had to assume the absence of magnetic 
monopoles and magnetic currents. Hence magnetic helicity is conserved in a highly conducting plasma provided there 
are no magnetic monopoles. Then we might expect that if magnetic helicity is externally driven to large values in 
a plasma of extremely high electrical conductivity, magnetic monopoles might be produced to cause the helicity to 
dissipate. 

Our comments above also apply to the creation of sphalerons in the electroweak model since the sphaleron is precisely 
a bound state of an electroweak monopole and antimonopole. The monopole pair in the sphaleron are confined by a 
Z-string [21] and cannot escape. Hence, if a sphaleron is produced, it will later decay. The production and decay of a 
sphaleron can lead to baryon number violation. So the discussion above suggests that it may be possible to observe 
baryon number violation by focusing magnetic helicity in a dense plasma 4 . 



3 In the SU(2) limit of the electroweak model there are no photons in the model. In this limit, sphaleron decay yields 8 Higgs particles 
and 42 W particles [20]. 

4 As a first step in further investigations it may help to study sphaleron production in the 1+1 dimensional toy model proposed in 
Ref. [22]. 
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While we are still speculating, it is possible that the magnetic field of certain atomic nuclei contain non-trivial 
magnetic helicity: since magnetic helicity violates parity, the lightest stable candidate nuclei are Li 7 and Be 9 [23]. 
The scattering of a large number of such nuclei, ~ 10 3 based on energy requirements, seems like a possible path to 
concentrate magnetic helicity in a small spatial volume, perhaps sufficient to form a sphaleron. One way to investigate 
this process theoretically would be to use the description of nuclei in terms of Skyrmions. To study the feasibility of 
baryon number violation in the collision of "helical nuclei" , we would need to couple the Skyrmion fields to electroweak 
gauge fields, then consider the scattering of corresponding "helical Skyrmions" in an attempt to find in-states that 
lead to sphaleron production, and thus baryon number violation. 

To conclude, we have studied the creation of solitons in a model that contains a twist degree of freedom in addition 
to the ones that are used for building the soliton. In our model in Sec. I, twist is an essential ingredient for the 
production of solitons because the dynamics in the zero twist sector is completely integrable. We have explicitly 
shown a range of parameters that lead to soliton creation. We have also discussed some applications of these ideas to 
the production of vortices and monopoles, and baryon number violation via the production of electroweak sphalerons. 
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Appendix A: Numerical Notes 

The most straightforward approach to solving the equations of motion is to use spherical coordinates as in Eqs. (9) , 
(10). The issue there is that spherical coordinates are not globally well-defined and one must use two different patches, 
say, (9 Z , 4> z ) with respect to the z— axis and (9 X , 4> x ) with respect to the a;— axis. If Z becomes small, we can solve the 
evolution equation in terms of X and <j) x , and then transform to 6 Z and 4> z . We found that this scheme introduced 
unacceptable numerical errors in the boundary region of the two patches and we decided not to pursue this scheme. 

Instead we decided to work with the vector representation (n) where the equation of motion can be written in first 
order form as 

9 t n = to (Al) 
dtfn = d x n — [to 2 — (d x n) 2 ]n — 713(71371 — e$) (A2) 

We evolved these equations using the explicit second order Crank-Nicholson method with two iterations [24] . We also 
implemented boundary conditions, though in practice our lattice was large enough that the boundaries were irrelevant. 

The issue now is that when we discretize the above equations, they do not preserve the condition n 2 = 1 and the 
orthogonality condition n ■ fn = 0. To remedy this, we enforced these constraints at every time step by rescaling n so 
as to get a unit vector, and by subtracting out from to the component parallel to h. To our satisfaction, this procedure 
conserved energy at the percent level and the evolution also remained stable. 
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